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Abstract. In this paper, we will analyse the cosmological models in Einstein-aether gravity, 
which is a modified theory of gravity in which a time-like vector field breaks the Lorentz 
symmetry. We will use this formalism to analyse different cosmological models with different 
behavior of the scale factor. In this analysis, we will use a certain functional dependence of 
the dark energy on the Hubble parameter. It will be demonstrated that the aether vector 
field has a non-trivial effect on these cosmological models. We will also perform the Om 
diagnostic in Einstein-aether gravity. Thus, we will fit parameters of the cosmological models 
using recent observational data. 
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1 Introduction 

The Lorentz symmetry is one of the most important symmetries in nature, and all 
particle physics experiments have demonstrated that this symmetry is not broken at the scale 
at which such experiments are performed. However, it is predicted from quantum gravity that 
the Lorentz symmetry should break down at Planck scale, where even the manifold structure 
of spacetime breaks down due to quantum fluctuations. In fact, almost all approaches to 
quantum gravity predict that the local Lorentz symmetry of spacetime only exists in some 
IR limit of the theory. So, the Lorentz symmetry is expected to break in the UV limit. It may 
be noted that it has been explicitly demonstrated that such a breaking of Lorentz symmetry 
in the UV limit occur in the discrete spacetime [1], models based on string field theory 
[2], spacetime foam [3], spin-network in loop quantum gravity (LQG) [4], non-commutative 
geometry [5, 6], and ghost condensation in perturbative quantum gravity [7]. As the Lorentz 
symmetry fixes the form of the energy-momentum dispersion relation, the breaking of Lorentz 
symmetry in the UV limit, will also lead to a modification of the energy-momentum dispersion 
relation in the UV limit. In fact, there are indications from the Greisen-Zatsepin-Kuzmin limit 
(GZK limit) t the usual energy-momentum relation will get modified in the UV limit [8, 9]. 
The Pierre Auger Collaboration and the High Resolution Fly’s Eye (HiRes) experiment have 
confirmed earlier results of the GZK cutoff [10]. So, it is possible that the Lorentz symmetry 
will break in the UV limit, and only occur as an effective symmetry in the IR limit. Thus, 
it is important to construct a theory, such that it will reproduce the general relativity in the 
IR limit, and break the local Lorentz symmetry in the UV limit. Such a theory has been 
constructed by using different Lifshitz scaling for space and time, and this theory is called 
the Horava-Lifshitz gravity [11, 12], This original proposal for the Horava-Lifshitz gravity 
improves the renormalization of gravity, as it differs from general relativity in the UV limit. 
However, there are several problems associated with this proposal, and these include the 
problems associated with instabilities, overconstrained evolution, and strong coupling at low 
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energies [13]-[18]. These problems occur due to a badly behaved scalar mode of gravity, which 
is produced by the presence of a nondynamical spatial foliation in the action. To resolve this 
problem, an extension of Horava-Lifshitz gravity called the BPSH theory has been proposed 
[19]. It has been demonstrated that the BPSH theory is equivalent to general relativity 
coupled to a dynamical unit timelike vector field [20]. Here the vector is restricted in the 
action to be hypersurface orthogonal. 

The phenomenology and observational constraints on the coupling parameters of Einstein- 
aether gravity have been studied [21], It may also be noted that constraints on Einstein-aether 
gravity from binary pulsars have also been discussed [22]. In this work, the consequences of 
Lorentz symmetry, which occur in Einstein-aether gravity, on the orbital evolution of binary 
pulsars. In the focus of this study was on the dissipative effects in such a process. It was 
observed that the breaking of the Lorentz symmetry modified such effects. Thus, the orbital 
dynamics of binary pulsars was also modified in Einstein-aether gravity. Such a modification 
causes the emission of dipolar radiation, and this made the orbital separation decrease faster 
than in general relativity. The quadruple component of the emission was also modified. The 
orbital evolution depends critically on the sensitivities of the stars, as this measure how their 
binding energies depend on the motion relative to the preferred frame. In this study such 
sensitivities have also been numerically calculated in Einstein-aether gravity. These predic¬ 
tions have been compared with observations, and this has been used to set constraints on 
Einstein-aether gravity. 

It has been demonstrated that the Einstein-aether theory can be analysed in the frame¬ 
work of the metric-affine gravity [23] . Such a formalism resembles the gauge theory of theory. 
In this formalism, the aether vector field is related to certain post-Riemannian nonmetricity 
pieces contained in an independent linear connection of spacetime. Black hole solution have 
also been studied in the Einstein-aether gravity. It has been demonstrated that the deviations 
from the Schwarzschild metric are typically only a few percent for most of the explored pa¬ 
rameter regions, and this makes it difficult to observe with electromagnetic probes, but they 
can be detected using gravitational wave detectors [24], As gravitational wave detectors are 
going to be used extensively in future astronomy, it is interesting to study the implications 
of Einstein-aether gravity. 

As the Einstein-aether gravity, introduce a time like vector field, they are expected to 
modify the cosmological evolution of the universe. In fact, various different solution for the 
accelerating universe in the Einstein-aether gravity have been studied [25]. These solutions 
have been used to analyse the inflationary behaviour of the early universe and late-time 
cosmological acceleration. It has been demonstrated that the aether field produces accelerated 
expansion in situations where inflation would not occur in general relativity. Hence, the 
aether held can effect the inflation in a very non-trivial way. The cosmological evolution of 
cosmological models based on Einstein-aether gravity with power-law potential have also been 
studied [26]. The cosmological models have also been studied in the Einstein-aether gravity 
coupled to a Galileon type scalar held [27]. It was observed that in such models, the universe 
experiences a late time acceleration, for pressure-less baryonic matter. It may be noted that 
gravitational wave can be used to analyse the cosmological aspects of Einstein-aether gravity. 
This is because it has been demonstrated that for cosmological models based Einstein-aether 
gravity, a direct correspondence exists between perfect fluids carrying anisotropic stress and 
a modification in the propagation of gravitational waves [28]. As the anisotropic stress can 
be measured in a model-independent manner, the gravitational waves can be used to obtain 
constraints on the cosmological models in the Einstein-aether gravity. Even though several 
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studies have been done on the Einstein-aether gravity, it is important to perform an extensive 
study of how this modification of general relativity can change cosmological models, and 
how it fits with the current data. This is important as many aspects of Einstein-aether 
gravity can be detected using gravitational waves, and in near future, it is expected that the 
gravitational wave will be used to study many of these interesting phenomena. So, in this 
paper, we perform a detail study on the modification of different cosmological models from 
Einstein-aether gravity. These models have been studied in general relativity, and we will, 
and we will analyse them in the Einstein-aether gravity. So, in this paper, we will first use 
the reconstruction technique to find some viable forms for Einstein-aether gravity. We obtain 
the expressions of the modified Friedmann equations and from these equations, we can find 
the effective density and the effective pressure for the Einstein-aether gravity. We will also 
fit the model with observational data. This will be done by using the cosmographic analysis 
involving the Om parametrization. We will also use the Snla, BAO and Hubble data to find 
the lcr and 2a contours for density parameter £l m arising from the Sne la + BAO. 


2 Einstein-aether Gravity 


In this section, we will review the main features of the Einstein-aether gravity. It may be 
noted that Einstein-aether gravity is equivalent to the BPSH generalization to the Horava- 
Lifshitz gravity [20], and so it breaks the Lorentz symmetry. However, the cosmology de¬ 
scribed by this theory would still be described by the standard Friedmann equations with 
an additional matter contribution [25]. This is because the breaking of Lorentz symmetry 
occurs due to a time-like vector field in the Einstein-aether theory, and so the cosmological 
effects can be obtained by analysing the correction to the standard Friedmann equations from 
this additional time-like vector field. The action S of the Einstein-aether gravity is given by 
[29, 30], 


S 



R 

~4nG 


+ Lea + L m j 


( 2 . 1 ) 


where Lea indicates the Lagrangian density of the aether vector field, while L m indicates the 
Lagrangian density of the usual matter fields. The Lagrangian density of the aether vector 
field Lea can be expressed as 


Lea = ^g f{k) + T^g x{a ° a -‘ + 1) ’ (2 - 2 > 

K = M~ 2 K^V a A c X7 b A d , (2.3) 

Kf d = c\g ab g cd + c 2 d a c S b + c 3 S a d S b , a, b = 0,1,2,3 (2.4) 


where ci, C 2 and c 3 are three dimensionless constant parameters, M is a coupling constant 
parameter, A is a Lagrangian multiplier, and A a is a contravariant vector. Here F ( K ) is an 
arbitrary function of the parameter K , and a function of the Hubble parameter H. 

Now using Eq. (2.1), the field equations, for this theory, can be written as [29, 30] 

q b rj-iEinstein—aether _|_ g 'jiQ r p rn (2 5) 

V a (F'J a b ) = 2A A b , (2.6) 


where 


F' 


T a — 
J h — 


dF 

dK’ 

2K a b d c V d A c . 


(2.7) 

( 2 . 8 ) 
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Here is the energy-momentum tensors for the matter field, and aether i s the 

energy-momentum tensor for the aether vector field, 


rpm 

ab 

h Einstein—aether 


= (p + p) u a u b + pg ab , 




J d A h — J„A b — J, 


(' ab R 


F' 


—Y(ab)F' + -g ab M 2 F + A A a A b . 


(2.9) 


( 2 . 10 ) 


Here p is the energy density and p is the pressure of the matter. Also u a is defined as 
u a = (1,0,0,0), and so it represents the four-velocity vector of the fluid, A a = (1,0, 0,0). It 
may be noted that it is represented by a time-like unitary vector. Now we define Y ab as 


Y ab — c i 


(V d A a ) (v<%) - (’ V a A d ) (V a A d ) 


( 2 . 11 ) 


The subscript (ab) indicates a symmetry with respect to the two indices. 

We consider a Friedmann-Robertson-Walker (FRW) metric, which is given by 


ds 2 = —dt 2 + a 2 (t) 


dr 2 


1 — kr 2 


+ r 2 (dd 2 + sin 2 9d(p 2 ) 


( 2 . 12 ) 


where a ( t ) represents the scale factor (which gives useful information about the expansion 
rate of the universe), t is the cosmic time, and k is the curvature parameter. Here its value 
can be —1, 0 or +1 corresponding to an open, a flat, or a closed Universe. The range of 6 
and 4> are 0 < 6 < it and 0 < (f> < 2ir. The four coordinates (r, t, 9, (j)) are also known as 
co-moving coordinates. 

Using Eqs. (2.3) and (2.4), we can easily obtain the following general expression for the 
parameter K [29, 30], 


K = 


3 eH 2 

~W r ' 


(2.13) 


where e is a constant parameter while H = a/a is the Hubble parameter. So, using Eq. (2.5), 
we obtain the Friedmann equations modified by the Einstein-aether gravity, 


e 




e Jt {Hpr) + {~ 2H + §) 


The conservation equation can now be written as 


8v tG 
— P ’ 

8irG (p + p). 


(2.14) 

(2.15) 


p + 2>H (p + p) = 0, 


(2.16) 


where the p indicates a temporal derivative of p. 

We denote the effective energy density in Einstein-aether gravity by pea, and the ef¬ 
fective pressure in Einstein-aether gravity by pea ■ So, we can rewrite Eqs. (2.14) and (2.15) 
as 




2 H + 


2k 


8v tG 1 

^ 3 Pea > 

8nG (p + p) + (pea + Pea) ■ 


(2.17) 


- 4 - 


(2.18) 












Therefore, comparing Eqs. (2.14) and (2.15) with Eqs. (2.17) and (2.18), we can write 


PEA = 3 SH* (r - U) . 
PEA = (f' - X) 


= -PEA - 


Pea 

3 H ' 


e [HF' + HF'^j 


Using Eq. (2.19), we obtain 


(2.19) 


( 2 . 20 ) 


F' -— 

2 K 


Pea 
3 eH 2 ' 


( 2 . 21 ) 


This is equivalent to the following master equation (using the expression of 3 eF[ 2 derived 
from Eq. (2.13)), 



Pea 
KM 2 ' 


( 2 . 22 ) 


Here that a prime ' indicates a derivative with respect to K , and so we have F' = 4^. 

Using the expressions for pea and given by Eqs. (2.19) and (2.20), we obtain that the 
EoS parameter u>ea for the Einstein-aether gravity, 


uea 


Pea 

Pea 


( HF' + HF'^j 

3 # 2 (F' ~ 5 )' 


(2.23) 


3 Models for Dark Energy 


The astrophysical data obtained from distant la supernovae, large scale structure, 
baryon acoustic oscillations, weak lensing and cosmic microwave background indicate the 
existence of dark energy [31-40]. In this paper, we will analyse the dark energy models using 
Einstein-aether gravity. The effective density and the effective pressure produced by Einstein- 
aether gravity can be used to generate dark energy, if the condition pea+^Pea < 0 is satisfied 
(i.e., if the strong energy condition is violated). Thus, we obtain 

2H 2 (f'~ > - (HF' + HF'^j . (3.1) 


This is the equation can be used to analyse dark energy in Einstein-aether gravity. So, to 
analyse the effect of dark energy on cosmological models in the Einstein-aether gravity, we 
can use F (K) and w Einstein-aether- It may be noted that the modified effective Friedmann 
equation can be written as 


87T Gpd ar ken e rgy ^-'Snj=n-^-( j ^) E, ...) 


d n f(X,Y .,-) 
dX ni dY n2 ... ‘ 


(3.2) 


where f(X,Y,...) is the matter part K action. For a model of dark energy , for example, 
the holographic dark energy, we have p d arkenergy = Pdarkener gy {H , H ,..). We can also write it 
as pdarkenergy = Pdarkenergy{X , Y,...). So, if we can solve the partial differential equations for 
/ (X, Y, ...), we can obtain the effect of dark energy on such cosmological models. 
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The generalized Nojiri-Odintsov Holographic dark energy models can be used to analyse 
the dependence of the dark energy on the Hubble parameter [59]. In fact, using the Granda- 
Oliveros model [41], which is a specific kind of Nojiri-Odintsov Holographic dark energy 
model, we can write, 


L go = ( aH 2 + PH} \ (3.3) 

where H = a/a is the Hubble parameter, and H is the temporal derivative of H. This model 
is characterized by two constant parameters a and (3. As the dark energy dominates the 
present cosmological epoch, and its contribution to cosmological epoch near the Big Bang 
was negligible (i.e. the amount of dark energy increased with the expansion of the universe), 
the energy density can be assumed to be a function of the Hubble parameter H and its 
temporal derivative. Such a dependence is characterized by these two parameters. There are 
other physical reasons which motivate the Granda-Oliveros model. This is because if the IR 
cut-off chosen is given by the particle horizon, the Holographic dark energy models are not 
able to produce an accelerated expansion for the present cosmological epoch. However, if 
the future event horizon is used an the IR cut-off, then the Holographic dark energy models 
have a problem with causality. It is possible to resolve both these problems by using Granda- 
Oliveros model [41]. It may be noted that in the limiting case {«,/?} = {2,1}, Lgo becomes 
proportional to the Ricci scalar curvature Lgo °c R- However, based on the observational 
data, such a choice is not physical. This is because the values for these parameters which 
best fit the observational data have been obtained [42], These value for a non-flat universe 
are given by [42], 

a = 0.8824^^80(1^+0.2213^^ p = 0 .5016t°;°^(la)!K(2a). (3.4) 

These value for a flat universe are given by [42], 

« = 0.8502iS(M t 0 0 f 0 ll(2a), (3 = 0.4 81 7™ 2 (la) l™(2a). (3.5) 

It may be noted that for Granda-Oliveros model, energy density po can be written as 

p D = 3c 2 (aH 2 + f3H^J . (3.6) 

In this paper, we will analyse the effect of Einstein-aether gravity on the cosmological evo¬ 
lution using the Granda-Oliveros model. We will analyse this for different models for the 
evolution of the scale factor of the universe, and analyse such a model for the observationally 
motivated values of a and f3. 

It may be noted the Granda-Oliveros model has been generalized to a Chen-Jing model 
[43]. In this cosmological model, the dark energy is a function of the Hubble parameter 
squared (i.e. H 2 ) and its first and second temporal derivatives of the Hubble parameter H 
and H , 


Pd = 3c 2 


a (^j +/3H + ~fH 2 , 


(3.7) 


where a, (3 and 7 represent three arbitrary dimensionless parameters. The inverse of the 
Hubble parameter, i.e. H , is introduced in the first of the three terms of Eq. (3.7), so that 
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each of these three terms have the dimensions. 

It may be noted that in the limiting case corresponding to a = 0, we recover the energy density 
of dark energy given by the Granda-Oliver os model [44], [45]. Furthermore, in the limiting 
case, a = 0 , = 1 and 7 = 2 , we obtain the expression of the energy density of dark energy 

with the IR cut-off proportional to the average radius of the Ricci scalar curvature, L oc R ^ 1 / 2 
(when k = 0). In this paper we will analyse this model for various cosmological model, with 
different evolution of the scale factor. We will obtain general expression for various parameters 
for this dark energy model, and they can be compared to various observational data. 

4 Cosmological Models 

In this section, we will analyse the behavior of various cosmological models in Einstein- 
aether gravity. These will correspond to different evolution of the scale factor. We will analyse 
them for the Granda-Oliveros model, using the values obtained from observation. We will 
also analyse them for the Chen-Jing model. 

Now we start by considering the the power-law cosmology. This cosmological model 
is interesting proposal for the evolution of the scalar factor, and it has been motivate by 
the existence of the flatness and horizon problems in the standard cosmology [46]. In this 
cosmological model, it is possible to assume the following form for the evolution of the scale 
factor, 


a (t) = aot m . 


(4.1) 


where ao is the present day value of a(t). It is also important to only consider m > 0, and 
this will produce an accelerating universe [46]. It may be noted that for m > 1, the power-law 
cosmology can solve the horizon problem, the flatness problem, and the problem associated 
with age of the early universe [47, 48]. The power-law cosmology has been used for analysing 
the cosmological behavior in modified theories of gravity [49, 50]. Now for Einstein-aether 
gravity, it is possible to analyse the Granda-Oliveros cut-off for power-law cosmological model. 
Thus, for a non-flat universe, we obtain, 


F{K) = 

For a flat universe, we obtain, 

F{K) = 


2c 2 i7(0.8824m - 0.5016) 


em 


2c 2 i7(0.8502m - 0.4817) 


em 


+ c x Vk. 


+ C { \/k. 


(4.2) 


(4.3) 


It is also possible to analyse the power-law cosmology using the Chen-Jing model, we have 

2 c 2 K [2a + n(n/3 + 7 )] 


F(K) = 


err 


+ c a Vk. 


(4.4) 


In this section, C\ will denote integration constants, for example, here C \, C 4 denotes the 
integration constants. 

It is also possible to consider a different kind of power law [51, 52], 


a(t) = ao(t s - t ) n , 


(4.5) 
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where n > 0 and t < t s . Such models have a future singularity at finite time, and this is 
denoted by t s . The Granda-Oliveros cut-off for such model can be analysed for both a flat 
universe and a non-flat universe. For a non-flat universe, we obtain 


F 2^(0.8 8 24 n + 0.5016) + 

era 


For a flat universe, we obtain 


F (K) = ^*(0-8502n +0.4817) + ^ 
era 


Now using Chen-Jing model, we obtain 


F(K) 


2c^[ 2« + "W + 7)] +C4 ^ 

era 2 


(4.6) 


(4.7) 


(4.8) 


It is also possible to analyse intermediate inflation in Einstein-aether gravity. The in¬ 
termediate inflation has been used to obtain exact analytic solutions for a given class of 
potentials for the inflation. The scale factor for intermediate inflation can be expressed as 
[53, 54], 


a(t) 



(4.9) 


where B > 0 and 0 < 0 < 1. For a non-flat universe, we have 


( KM 2 \ 2 (-i+e) 

F(A, = W W(^) 


X 


o ( KM 2 \ 2 (-i+ 0 ) 

-0.5016(—1 + Of + 0.882450 f g2 j 


+Vkc 7 . 


(4.10) 


For a flat universe, we have 


i n ( KM 2 \ 2 (-i+») 
F( K ) = 2(^,-Va( 51w ) 

—0.4817(—1 + 6) 2 + 0.85020 


/ KM 2 \ 2 (-!+ e > 
\35 2 e0 2 y 


+Vkc 7 . 


(4.11) 

















Now for Chen-Jing model, we obtain 


, as( i w 2 ® oi-ef i \ 9 n ( KM 2 \ i 
F(K) = £ _1 2 3 4 Ui+eJ 2 i-i +9 J C 2 A -/ j 

r,0o( 1 \ 29 ( KM 2 \ 2i 

2 p V^r+ey VbwJ 

+ + *(=&»)” 

2 » 3 2-*(^F,)b('a_r 1+2l ’('iP£i 2 " ,( ^ ) ' 


+- 


+- 


T l -1+0 ) \b 2 e6 2 J 


2 1 + 3e B9 + 2 e B 


i 


- 1+0 


2 e 


9 + 2B 


l 


-i+0 


30 


\-=t+e) e a ( 2 - 39 + 9 2 ) 


B 2 


1 + 2 !+ 2e 


i 


- 1+0 


-20 


9 2 


, 


+c % Vk. 


(4.12) 

(4.13) 


It is possible to analyse universes with no past time-like singularity. In such cosmological 
models, the universe in the infinite past there exists a static state of cosmology, and this 
state evolves into an inflationary stage. The scale factor for such cosmological models can be 
expressed as [55, 56], 

a(t) = A(B + e nt ) X (4.14) 

where A, B , n and A are four positive constant parameters. In order to avoid singularities, 
we have to use B > 0. Furthermore, for for the positivity of scale factor, we have to use 
A > 0. It may be noted that in this model, for a < 0 or A < 0, a singularity exists. So, for 
the expanding model, we have to only consider a > 0 and A > 0. Using the Granda-Oliver os 
cut-off for this model, we can analyse a flat universe and a non-flat universe. So, for a non-flat 
universe, we have 


c 2 A^ M.7648A - 1.0032 + 0.5016y^g|^ log if] 

F(A 0 = --£- Z 

+VKC 5 . (4.15) 

For a flat universe, we have 

c 2 K f 1.7004A - 0.9634 + 0.4817y^JpTog AT J 
A W = --^--- Z 

+VKC 5 . (4.16) 

We can also use the Chen-Jing model, and obtain 

F (AT) = c 2 (eM 2 A 2 ) -1 [-6en 2 aA 2 + 2KM 2 (2a - 7 A + /3A 2 ) 

+ V3KM 2 ^^(-3a + 1 \)logK , 

+c 6 Vk. 
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It is possible to analyse matter dominated universe and the accelerated phase of the 
universe using a single formalism [57, 58]. In such cosmological models, the Hubble constant 
is given by [59, 60] 

H(t) = H 0 + ^-, (4.18) 

where Hq and Hi are two constant parameters. In this case, for the non-flat universe, we 
obtain 


F (K) = -2c 2 (9£M 2 H 1 )~ 1 [-9m88£ 2 KM 2 H 0 - 13.5431e.ffg 
+KM 2 (0.5016e 3 KM 2 - 7.9416ffi) 

+y/KC w . (4.19) 

Furthermore, for the flat universe, we obtain 

F{K) = -2c 2 (9£M 2 H 1 )- 1 [-8.6706e 2 KM 2 H 0 - 13.0059effg 
+KM 2 (0A817£ 3 KM 2 - 7.5618Hi)] 

+VKC W . (4.20) 


Now using the Chen-Jing model, we obtain 

2c 2 M 2 a K 3 / 2 c 2 M 2 {3K 3 / 2 


F(K) = VK 


C n + 


eHi 2 eH x 


( c 2 M 2r ) K 3 / 2 9c 2 V2MH 0 aK | 3c 2 V2MH 0 (3K 

^ ■) “ 1 “ ” 


+ 


e \/7 H i “ II i 

18^2 ec 2 H 0 2 Ha In (K) 3V2ec 2 H 0 2 H/3 In (. K ) 


H\ 2 M 


H\M 


6\/2ec 2 a Hq 3 In (K) 
H x 2 M 


(4.21) 


The quantum deformed de Sitter (g-de Sitter) solution has been obtained by a quantum 
deformation of the quantum deformation of the conformal group [61]. In fact, the g-deformed 
de Sitter solution has also been used in the analyzing of dS/CFT correspondence and the 
entanglement entropy for such a solution has also been obtained [61]. The g-de Sitter has also 
been used for analyzing cosmology [62]. Now we will analyse this model for the g-de Sitter 
scale factor [62], 


a(t) = e q (H 0 t) 


1 + (g - l)H 0 t 


1 

q~ 1 


(4.22) 


In this model it is possible to interpolate between the cosmological model based on a power- 
law and the cosmological model involving de Sitter spacetime. In fact, for early times, Hot S> 
1, we obtain 


Nearly (t') ~ 



t p . 


(4.23) 
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It is possible to have an acceleration expansion, when p > 1 , and g < 2 . So, we can write 

a e arly(t ) + e q (H 0 t ) + a d 5 (t). (4.24) 

This inequality, given in Eq. (4.24), produces interesting cosmological evolution in g-de Sitter. 
The g-de Sitter can be used to smoothly connect the early cosmological epoch to late time 
evolution universe. Now we can analyse such a model for a non-flat universe as 


F(A') = 6 c 2 Hi 1 K ' 




2 +- 


9 (- 2 +q) log[K]+log 


£ H 2 

£n 0 


—Y 


+ 2 ( — 2+q) log 


1 + 9 


-1 + 9 


2( 2+9)+ 1+9) 


— 1 + 9 

K \ 2 ( —2 +q) 

^ 0 


xK 4-69+2 9 2 

(-2+ g ) 2 -0.5016 


A 


^0 


x — 


-1 + 9 

We can also analyse such a model for the flat universe as 


-1+29 

2( — 2+9( — 1 +9) 

0.8824 

M~ 2 + \[KC\\. 


F(K) = 6 c 2 H 2 


K \ 1+ ^+7 


eH 2 ) 


eH 2 


K 


— 1+9 


9 (- 2 + 9 ) log[fC]+log 


K 


\ — ~^-+9 
K \ 2 ( —2-fq) 

7h$ 


XK 4 - 69 + 2 < j 2 

(-2 + q) 2 -0.4817 


K 


eH 2 


x — 


A + 2 (— 2 + 9 ) log 
eH l 

2 ( —2 + 9 ) 2(—1 + q ) 


-1 + 29 

2 (- 2 + 9 )^(-l+ 9 ) 

0.8502 

M -2 + VkCh. 


-1 + 9 

We can use the Chen-Jing model, and obtain 

, . ,— 3c 2 (2 a g 2 — 4 a g + 2 a — (3 q + f3 + 7 ) 

F(A) = ^KC 13 + -1-----—— - 


(4.25) 


(4.26) 


(4.27) 


Thus, we can see the Einstein-aether gravity modifies the cosmological evolution in various 
different model, with different evolution of scale factor. Hence, this deformation is almost 
a universal feature of Einstein-aether gravity. We have calculated F(K ) for these different 
cosmological models. Thus, these cosmological model directly depend on the aether vector 
field. It is also possible to calculate Lqo , Pea and pea for these various different 

cosmological models. It can be argued that these quantities will also depend on the aether 
field (see Appendix). Hence, the breaking of Lorentz symmetry by the introduction of a time¬ 
like aether vector field can modify the cosmological dynamics in a non-trivial way. Here we 
explicitly calculated such a modification for a large number of cosmological models. 
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5 Om Diagnostic Analysis 


In this section, we will perform the Om diagnostic analysis of various different cos¬ 
mological models. The cosmological parameters like the Hubble parameter H, deceleration 
parameter q, and the Equation of State (EoS) parameter oj are important to understand the 
behavior of cosmological models. It is theoretically and observationally known that different 
dark energy models produce a positive Hubble parameter and a negative deceleration pa¬ 
rameter (i.e. H > 0 and q < 0, for the present cosmological epoch. So, H and q can not be 
used to effectively differentiate between the different dark energy models. Therefore, a higher 
order of time derivatives of a(t) is required to analyse the dark energy models [65, 66]. So, 
third order temporal derivative of a(t) can be used to resolve the problem that most dark 
energy models produce H > 0 and q < 0 for the present cosmological epoch. So, now the 
statefinder parameters {r, s}, can be expressed as 

(5.1) 

(5.2) 

where q represents the deceleration parameter, which is given by 


r = 


s = 


r — 1 

3 (<7-1/2)’ 


1 a 

An alternative way to write r and s is as 


r 

s 


l + 2, m + m' 

3 HH + H 
3 H (2H + 3H* 


3H + H/H 
3 ( 2H + 3f7 2 ) 


(5.3) 


(5.4) 


(5.5) 


It may be noted that statefinder parameters {r, s} = {1,0} represents the point where the 
flat ACDM model exists in the r — s plane [67]. So, the departure of dark energy models 
from this fixed point can be used to obtain the distance of these models from the flat ACDM 
model, taken as reference model. 

We also note that in the {r, s} plane, a positive value of the parameter s (i.e. s > 0) indicates 
a quintessence-like model of dark energy, and a negative value of the parameter s (i.e. s < 0) 
indicates a phantom-like model of dark energy . Furthermore, the evolution from phantom 
to quintessence is obtained by crossing of the point {r, s} = {1,0} in the {r, s} plane [68]. 
So, different cosmological models, like the models with a cosmological constant A, braneworld 
models , chaplygin gas and quintessence models, have been studied using such an analysis 
[66]. In this study, it was argued that {r, s} can be used to differentiate between different 
models. An analysis based on {r, s} has also been used to differentiate between dark energy 
and modified gravity [68, 69]. 

An important geometrical diagnostic which can be used to for such analysis is called the 
Om diagnostic analysis [70]. Usually in the study of the statefinder parameters r and s, 
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Figure 1. r — s, r — q, m — iv Einstein-aether for power-law in the redshift t = (1 + z) 1//m range 
0.07 < z < 2.3. 


higher order temporal derivatives of a(t) are used. However, in the Om diagnostic analysis 
only first order temporal derivative are used. This is because it only involves the Hubble 
parameter, and the Hubble parameter depends on a single time derivative of a(t). So, the 
Om diagnosis can be considered as a simpler diagnostic than the statefinder diagnosis [71]. 
It may be noted that the Om diagnosis has also been applied to Galileons models [72, 73]. 
This set of parameters can now be represented as 


Om(z ) 


HQ) 1 

H 0 \ 1 

(1 + z ) 3 — 1 ' 


For a constant EoS parameter u, the expression for Om(z ) is given by 

(1 _|_ 2 ) 3 ( 1 + w ) _ i 

Om(z) — H m o + (1 fi m o) — ■ To ■ , 

(1 + z) 6 - 1 


(5.6) 


(5.7) 


Thus, we observe that we have different values of Om(z) = O m o for the ACDM model, 
quintessence, and phantom cosmological models. In Figure 1, we plot the first cosmological 
parameters r — s, r — q, m — oo Einstein-aether for power-law in the redshift t = (1 + z) _1 / m 
range 0.07 < z < 2.3. A continuous behavior is observed. For r — s, we observe that when r 
is increasing, s is starting to decreasing monotonically, never vanishes. A similar pattern is 
repeating but in the negative range for q. As we observe, the observational value for q ~ —0.67 
exists in this model. In Figure 2, we plot Om(z) for power-law in the redshift t = (1 + z) _1 / m 
range 0.07 < z < 2.3. We observe that when the redshift z is increasing within the interval 
0.07 < z < 2.3, the Om(z) is decreasing monotonically. For all power exponents m ^ 2, 
always Om(z) > Q m o, so the model mimics a quintessence with effective EoS w > —1. 

In Figure 3, we plot the first cosmological parameters r — s, r — q, oj Einstein-aether for 
a cosmological model with a future singularity. In the redshift t = t s — (1 + z) 1 /” range 


- 13 











z 


m=l.m=2-m=3 — ■ — m=4 


Figure 2. Om(z) for power-law in the redshift t = (1 + z) x ' m range 0.07 < z < 2.3. 


0.07 < * < 2.3, a continuous behavior is observed. For r — s, we observe that when r 
is increasing, s is starting to decreasing monotonically, always remains negative. A similar 
pattern is repeating but in the negative ranges for r, q. As we see, the observational value for 
q ~ —0.67 do not exist in this model. In Figure 4, we plot Om(z) for future singularities model 
in the redshift t = t s — (1 + z )■*■/" range 0.07 < z < 2.3. We observe that when the redshift z is 
increasing within the interval 0.07 < 2 < 2.3, the Om(z) is decreasing monotonically, always 
Om{z ) > fi m o- So, the model with future singularity mimics a quintessence with effective 


EoS w > —1. 


In Figure 5, we plot the first cosmological parameters r — s, r — q, uj Einstein-aether for 

In (A + zA) 

X ~ B 


n 


-l 


range 0.07 < z < 2.3. 


models of emergent universe, in the redshift t = In ^ 

A continuous behavior is observed. For r — s, we observe that when r is increasing, s starts 
to increase monotonically too, always remaining positive. A similar pattern is repeating for 
r, q. As we observe, the observational value for q ~ —0.67 do not exist in this model. In 


In (A + zA) 

X — 


n 


-l 


Figure 6 , we plot Om(z ) for emergent universe in the redshift t = In 
range 0.07 < z < 2.3. We observe that when the redshift 2 is increasing within the interval 
0.07 < 2 < 2.3, the Om(z) is monotonically decreasing, and Om(z) < Q m o, so this model 
mimics a phantom model with effective EoS w < — 1. 

In Figure 8, we plot the first cosmological parameters r — s, r — q, uj Einstein-aether for for 


intermediate inflation, in the redshift t = 


ln(l+z) 
B 


range 0.07 < z < 2.3. A continuous 


behavior is observed. For r — s, we observe that when r is increasing, s = — 1 remains 
constant. For r, q when we are increasing r, q is decreasing and —1 < q < 0. As we observe, 
the approved observational value for q ~ —0.67 do not exist in this model. In Figure 8 , we 


1 1 


plot Om(z) for intermediate inflation in the redshift t = — range 0.07 < 2 < 2.3. 

We observe that when the redshift 2 is increasing within the interval 0.07 < 2 < 2.3, the 
Om(z) is monotonically increasing, and Om(z) < Q m 0 , so the intermediate inflation mimics 
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r-s.r-q- n-(ss EA 


Figure 3. r-s, r-q, toEinstein-aether for model with future singularity in the redshift t = t s — (1 + z) 1 ^ 71 
range 0.07 < z < 2.3. 


Figure 4. 

x < 2.3. 


Om(z) 



■ n=l 


n=2-n=3-n=4 


Om(z ) for model with future singularity in the redshift t = t s — (1 + z) 1 / n range 0.07 < 


a phantom model with effective EoS w < — 1. 

6 Observational Constraints 

In this section, we will apply observational data from la Supernovae la, baryonic acous¬ 
tic oscillations (BAO), and data of the Hubble parameter H to study the constraints on 
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r-q 


( InfA+zA) _\ 

e * — B j n 

range 0.07 < z < 2.3 with parameters A = 1; B = 1; A = 2; n = 1. 



Z 


, , ( ln(A + zA) \ 1 

Figure 6. Om[z) for emergent universe in the redshift t = lnle ^ —Bjn range 0.07 < z < 
2.3 with parameters A = 1; B = 1; A = 2; n = 1. 


parameters of different cosmological models. The total x' 2 f° r joint data set which we use is 
defined by 


Xtot ~ Xsn + Xbao + XHub 


( 6 . 1 ) 


where the xf f° r each set of data is evaluated. To compute it we need the luminosity distance 

D L (z). 
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Figure 7. r-s, r-q, uiEinstein-aether for intermediate inflation in the redshift t = 
0.07 < z < 2.3 with parameters 0 = 2; B = 1. 


ln(l+z) 

B 


e - 1 

range 



Z 


Figure 8. Om(z ) for intermediate inflation in the redshift t 
with parameters 9 = 2; B = 1. 


ln(l+z) 

B 


e- 1 

range 0.07 < z < 2.3 


The luminosity distance is defined as 


Dl{z) 



Hodz' 

H&)' 


( 6 . 2 ) 
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Table 1. Values of ottzslo) l° r distinct values of Zbao [63]. 


ZBAO 

±l( 2 *) 

Dv(Zbao) 

0.106 

30.95 ± 1.46 

0.2 

17.55 ±0.60 

0.35 

10.11 ±0.37 

0.44 

8.44 ±0.67 

0.6 

6.69 ±0.33 

0.73 

5.45 ±0.31 


We use the distance modulus fi, which is given by 

M = m - M = 51og.Dk + no, (6.3) 


where m and M are defined as the apparent and absolute magnitudes of the Supernovae. 
Here fiQ = 5log + 25 is a nuisance parameter (which will be marginalized). We then 

have that the corresponding x 2 for this data set, 


x!n(Mo, 0 ) — 


^ ^ [l^th^Zii MOi ^) Mobs^Zi)] 
2=1 




(6.4) 


where fj, 0 b s , Mth and indicates the observed distance modulus, the theoretical distance mod¬ 
ulus and the uncertainty in the distance modulus, respectively. Furthermore, the parameters 
in the cosmological models are indicated by 9. For example, for the power law reconstruction 
scheme it given by m, the exponent of in the q -de Sitter it is given by the non-extensivity 
parameter q. Now we obtain, 


where 


XsnW — MO) 


B{9) 2 


580 

^w = E 

i— 1 
580 

m = E 
2=1 
580 

cm = 


[Mth (Zi > A^O 0;^) Mobs (-^i)] 

0 ±(Zj ) 2 

M0 — 0 ) 9) flobs^Zi') 
cr^{zi) 2 

1 

c y^Zi) 2 ' 


(6.5) 


( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 


If we use BAO data of Dvfz^Ao) ’ we ^ ave z * ~ 1091 as the decoupling time, (Ia(z) = fj u{z') 

1 

as the co-moving angular-diameter distance and Dy(z) = [dA{z) 2 j^-^j 3 as the dilation 
scale. Using this data set, the x'bao is defined as 


X 2 bao = X T C~ 1 X. 


(6.9) 
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Here what is needed is in the following column vector, 


X = 


( 


\ 


<4(z*) 

LV(0.106) 

D v ( 0.2) 
dA(z+) 

D v { 0.35) 
d-A (z*) 

D V (0A4) 

dA(z*) 

D v { 0.6) 
dAjz*) 

D v { 0.73) 


- 30.95 \ 
17.55 
10.11 
-8.44 
6.69 
-5.45 ) 


( 6 . 10 ) 


Furthermore, the C _1 is the inverse covariance matrix. Finally, we use the observational data 
on Hubble parameter as recently compiled by [63] in the redshift range 0.07 < z < 2.3. In 
this data set, the Hubble constant Hq is taken from the PLANCK 2013 results [64], It may 
be noted that the the normalized Hubble parameter is defined by h = H/Hq. In this data 
set, the x 2 f° r the normalized Hubble parameter is computed as 


29 

XHub(0) = 


i=l 


[hth(zi,6) h 0 bs(^j)] 
&h(zi ) 2 


( 6 . 11 ) 


where h D b s is the observed value of the normalized Hubble parameter, and is theoretical 
values of the normalized Hubble parameter. The error can now be estimated as 


<?h = 


<?H VHa 

H H 0 


h, 


( 6 . 12 ) 


where an is the error in H, and ajj 0 is the error in Hq. 



Figure 9. This figure shows the la (plotted in dark) and 2 a (plotted in light) likelihood contours for 
different cosmological with joint data (Snla+Hubble+BAO). 

In Figure 9, we plot the la (dark regions) and 2cr (light regions) likelihood contours 
for these cosmological models, Using the joint data (SNIa+Hubble+BAO), we observe that 
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the best fit value of the parameters which are found to be = 0.319 . Thus, for models 
with a power-law the best fit occurs for m = 3.218lQQ 5 g4(la) +qF urthermore, it 
is possible to have analyse certain models with a future singularity after finite time, and 
for these models, the best fit occurs for n = 4.017+° 0I53 (1°") -o 0876(2' 7 )- The best fit for 
emergent universe occurs for n = 2.054+°;°|^(1ct) A = 6^;^(la)+^||^(2cr), 

and the best fir for intermediate inflation occurs for B = 2 . 036 ^qq 2 ii( l* 7 )-o 0465 (^ a )’^ = 
°-756 ±o;o765 (la) ±oS 2(7 )- Thus, we have analysed different cosmological models in Einstein- 
aether gravity, and used observational data to analyze the value of parameters in these cos¬ 
mological models. 


7 Conclusions 

In this paper, we analysed various different cosmological models based on the Einstein- 
aether gravity. In Einstein-aether gravity, a time-like vector field couples the usual Einstein 
Lagrangian, and this time-like vector held breaks the Lorentz symmetry of the theory. In 
this paper, we have analysed various different cosmological models using Einstein-aether 
gravity. It was demonstrated that the aether held modihes the cosmology in a non-trivial 
way. Explicit expressions for such a modification to various different cosmological models 
were derived in this paper. Furthermore, the cosmological models based on Einstein-aether 
gravity were also compared with observational data. This was done by using the cosmographic 
analysis involving the Om parametrization. Thus, the Snla, BAO and Hubble data was used 
to obtain the la and 2cr contours for density parameter fl m arising from the Sne la + BAO. 

It is important to perform such an analysis as it is expected that gravitational waves 
can be used to test Einstein-aether gravity, and as gravitational wave will be used to test 
several of the predictions of Einstein-aether gravity, in near future, it is important to analyse 
the effect of Einstein-aether gravity on cosmology. In fact, it has been predicted that grav¬ 
itational wave detectors can be used to test Einstein-aether gravity [24]. Thus, it becomes 
important to analyse various different cosmological models using Einstein-aether gravity. It 
may be noted that as the Einstein-aether gravity modifies the cosmological models in a non¬ 
trivial way, it would also be interesting to analyse quantum cosmology using these modified 
cosmological models. It would be possible to calculate the Wheeler-DeWitt equation for these 
cosmological models, and the wave function of the universe can then be obtained as a so¬ 
lution to the Wheeler-DeWitt equation. We would like to mention, that such an analysis 
would be very interesting and important. Furthermore, as the time-like vector field breaks 
the time-reparametrization symmetry, it would modify the Wheeler-DeWitt equation in a 
very non-trivial way. It might be possible to use this time-like aether vector field to obtain 
a direction of time, even in the Wheeler-DeWitt equation. Thus, it might be possible that 
this formalism can be used as a solution to the problem of time. It would be interesting to 
perform such an analysis, for these cosmological models. 

It may be noted that the Horava-Lifshitz gravity has been used for analyzing type 
IIA string theory [75], type IIB string theory [76], AdS/CFT correspondence [77-80], dilaton 
black branes [81, 82], and dilaton black holes [83, 84]. As the Horava-Lifshitz gravity is related 
to the Einstein-aether gravity [20], it would be interesting to analyse these systems using 
Einstein-aether gravity. In fact, it has been demonstrated that Einstein-aether gravity can be 
related to the noncritical string [74]. Thus, it would be interesting to analyse this connection 
further, and also study various cosmological models motivated from string theory in Einstein- 
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aether gravity. It may be noted that the Einstein-aether gravity has been demonstrated to 
be equivalent to generalization of Horava-Lifshitz gravity 


8 Appendix 


In this appendix, we will explicitly calculate various cosmological solutions in Einstein- 
aether gravity. 

The first model which we are studying is the power law, 


a (t) = aot m . 


( 8 . 1 ) 


where ao is the present day value of a (t ), and we must have m > 0 for an accelerating 
universe [46]. With the choice of scale factor made in Eq. (8.1), we obtain that the Hubble 
parameter H, 


H = 


a (t) m 

a (t) t 


( 8 . 2 ) 


Moreover, we have that the first and the second time derivatives of the Hubble parameter 
obtained in Eq. (8.2), 


■ rri 

H t 2 ' 

(8.3) 

2 m, 

H = 

t 3 

(8.4) 


Furthermore, using in Eq. (2.13) the expression of H derived in Eq. (8.2), we obtain that the 
expression of K, 


K = 


3 em 2 
M 2 t 2 ' 


(8.5) 


Using in the general expression of Lqo given in Eq. (3.3) the expressions of H and H obtained 
in Eqs. (8.2) and (8.3), we have, 


Lgo = 


y/m(ma — (3) 


( 8 . 6 ) 


Therefore, we can conclude that the expression of pea with Granda-Oliveros cut-off can be 
written as 

3 c 2 m(ma — /3) 


Pea = 


t 2 


(8.7) 


Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.2) and (8.7) or equivalently 
in Eq. (2.22), the expressions of K and pea given in Eqs. (8.5) and (8.7), we obtain the 
following differential equation for F(K), 


dF (K) F(K) (?{ma — /3) 


dK 

which has the following solution, 


2 K 


= 0 , 


£171 


F(K) = 


2c 2 K (mot — /3 ) 


em 


+ CiVK, 


( 8 . 8 ) 


(8.9) 
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where C\ represents an integration constant. 

Using in Eq. (2.20) the expression of p derived in Eq. (8.7) along with the expression of H 
obtained in Eq. (8.2), we can write the pressure pea as 


Pea 


c 2 (2 — 3m) (ma — (3) 

T 2 


Therefore, the EoS parameter cue A is given by 


PEA 

uea =- 

Pea 


At Ricci scale, for a = 2 and (3 = 1, we obtain, 



Lgo 

Pea 
F (K) 

PEA 


t 

y/m(2m — 1) ’ 

3c 2 m(2m — 1) 

2Af(L— 1) ^ 

em 

<? (2 — 3m) (2m — 1) 

7 ' 


( 8 . 10 ) 

( 8 . 11 ) 

( 8 . 12 ) 

(8.13) 

(8.14) 

(8.15) 


Moreover, for a ~ 0.8824 and (3 ~ 0.5016, i.e., for the values of a and (3 corresponding to a 


non-flat Universe, we obtain, 



Lgo = 

t 

(8.16) 

v / m(0.8824m - 0.5016)' 

3c 2 m(0.8824m - 0.5016) 

Pea = 

(8.17) 

t 2 

F(K ) = 

2c 2 _fT(0.8824m - 0.5016) ^ c . ry 

(8.18) 

em 

Pea = 

c 2 (2 - 3m) (0.8824m - 0.5016) 

(8.19) 

t 2 


Furthermore, for a ~ 0.8502 and (3 ~ 0.4817, i.e., for the values of a and /3 corresponding to 
a flat Universe, we obtain, 


Lgo = 

t 

(8.20) 

y / m(0.8502m - 0.4817)' 

3c 2 m(0.8502m - 0.4817) 

Pea = 

(8.21) 

t 2 

F(K) = 

2c 2 7^(0.8502m - 0.4817) | ^ ^ 

(8.22) 

em 

Pea = 

c 2 (2 - 3m) (0.8502m - 0.4817) 

(8.23) 

t 2 


We now consider the Chen-Jing model studied in this paper, i.e., the one with the first 
and the second time derivatives of the Hubble parameter H. Using the expressions of H, H 
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and H given in Eqs. (8.2), (8.3) and (8.4) in Eq. (3.7), we obtain the expression of pea with 
higher derivatives of the Hubble parameter, 


Pea 


3c 2 [2a + m{m(3 — 7 )] 

T 2 


(8.24) 


Using the expressions of H and pea given in Eqs. (8.2) and (8.24) in Eq. (2.21) or equivalently 
in Eq. (2.22), the expressions of K and pea from Eqs. (8.5) and (8.24), we obtain the following 
differential equation for F (K), 


dF (K) F (K) c 2 [2a + m(m/3 - 7 )] 
dK 2 K em 2 

whose solution is given by 


F(K) 


2c 2 K [2a + m{mf3 - 7 )] | 

em 2 


(8.25) 


(8.26) 


where C 2 is an integration constant. Substituting in Eq. (2.20) the expression of pea derived 
in Eq. (8.24), along with the expression of H obtained in Eq. (8.2), we can write the pressure 
Pea as 


Pea = - 


c 2 (—2 + 3m) [2a + m {m(3 — 7 )] 
mt 2 


Therefore, the EoS parameter ujea is given by 


(8.27) 


PEA 

uea =- 

Pea 



(8.28) 


which is the same as Eq. (8.11). 

We now consider the second scale factor considered in this work, which is another form 
of the power law [51, 52], 


a{t) = a 0 {t s - t ) n , 


(8.29) 


where n > 0 and t < t s . Here a 0 is the present day value of a(t), while t s is the probable 
future singularity at finite time. So, this model has a future singularity. With the choice of 
scale factor made in Eq. (8.29), we obtain the Hubble parameter H, 


a ( t ) n 

a (t) t s —t 


(8.30) 


Moreover, the first and the second time derivatives of the Hubble parameter are given by 


H = 
H = 


n 

{t s ~t) 2 ' 
2 n 

{ts - t) 3 ' 


(8.31) 

(8.32) 


Furthermore, using in Eq. (2.13) the expression of H derived in Eq. (8.30), we obtain that 
the expression of K, 


K = 


Sen 2 

AP{t s -t) 2 ' 


(8.33) 
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Using in the general expression of Lgo given in Eq. (3.3), the expressions of H and H 
obtained in Eqs. (8.30) and (8.31), we obtain 


Lgo 


(t s ~ t) 

y/n(na + (3) 


(8.34) 


Therefore, we can conclude that the expression of pea with Granda-Oliveros cut-off can be 
written as 


Pea 


3c 2 n(na + f3) 

(t - t s ) 2 


(8.35) 


Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.30) and (8.35), or equiva¬ 
lently in Eq. (2.22) the expressions of K and pea given in Eqs. (8.33) and (8.35), we obtain 
the following differential equation for F (K), 


dF F c 2 (na + /3) 
dK 2 K en 


whose solution is given by: 


F(K) 


2 <?K{na + P) +c ^ 
en 


(8.36) 


(8.37) 


where C 3 represents an integration constant. 

Using in Eq. (2.20) the expression of pea derived in Eq. (8.35), along with the expression 
of H obtained in Eq. (8.35), we can write the pressure pea as 


Pea 


c 2 (2 + 3n) (not + f3) 

(t s ~ tf 


(8.38) 


Therefore, the EoS parameter ujea is given by 


uea 


Pea _ _2_ 

Pea 3 n 


At Ricci scale, i.e., for a = 2 and (3 = 1, we obtain 


Lgo 

Pea 

F(K) 

Pea 


(t s -1 ) 

y/n(2n + 1 ) ’ 

3c 2 n(2n + 1) 

(t - t s ) 2 

2 Af( 2 n + l) ^ 

en 

c 2 (2 + 3 n) (2n + 1 ) 

(t s ~ t) 2 


(8.39) 


(8.40) 

(8.41) 

(8.42) 

(8.43) 
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Moreover, for a ~ 0.8824 and j5 « 0.5016, i.e., for the value of a and j3 corresponding to a 
non-flat universe, we obtain 


Lgo 

Pea 
F(K ) 

PEA 


_ (t s ~ t) _ 

^/n(0.8824n + 0.5016) ’ 
3c 2 n(0.8824n + 0.5016) 

(t-t s ) 2 

2c 2 iv (0.8824n + 0.5016) | ^ ^ 
en 

c 2 (2 + 3n) (0.8824n + 0.5016) 

( t s - t ) 2 


(8.44) 

(8.45) 

(8.46) 

(8.47) 


Furthermore, for a ~ 0.8502 and j3 ~ 0.4817, i.e., for the value of a and j3 corresponding to 
a flat Universe, we obtain 


Lgo 

Pea 
F (K) 
Pea 


_ (t s - t) _ 

^/n(0.8502n + 0.4817) ’ 
3c 2 n(0.8502n + 0.4817) 

(t ~ t s ) 2 

2c 2 iv (0.8502n 4 - 0.4817) | ^ ^ 
en 

c 2 (2 + 3n) (0.8502n + 0.4817) 

(t s ~ tf ' 


(8.48) 

(8.49) 

(8.50) 

(8.51) 


We now consider the Chen-Jing model studied in this paper, i.e., the energy density 
model with the first and the second time derivatives of the Hubble parameter H. Using in 
Eq. (3.7) the expressions of H , H and H obtained in Eqs. (8.30), (8.31) and (8.32), we obtain 
the expression for pea, 

3c 2 

Pea = ^ _ t y [2a + n(nf3 + 7 )] . (8.52) 

Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.30) and (8.52), or equiva¬ 
lently in Eq. (2.22) the expressions of K and pea given in Eqs. (8.33) and (8.52), we obtain 
the following differential equation for F (K), 


dF F c 2 [2a + n(n/3 + 7 )] 
dK~ 2K en 2 


which solution is given by, 


(8.53) 


en z 

where C 4 represents an integration constant. 

Using in Eq. (2.20) the expression of pea derived in Eq. (8.52) along with the expression 
of H obtained in Eq. (8.30), we can write the pressure pea as follows, 


Pea = ~ 


c 2 (2 + 3 n) [ 2 a + n(nj3 + 7 )] 
n(t- t s ) 2 


(8.55) 
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Therefore, the EoS parameter ujea is given by, 


UEA 


Pea _ _2_ 

Pea 3 n ’ 


(8.56) 


which is the same results of Eq. (8.39). 

We can also analyse an emergent universe using this analysis. The scale factor for such 
a cosmological model is given by [55, 56], 


a(t) =A(B + e nt ) X 


(8.57) 


where A, B, n and A are four positive constant parameters. With the choice of scale factor 
given in Eq. (8.57), we can obtain the Hubble parameter H, 


H = 


e nt n\ 

B + e nt ' 


(8.58) 


Moreover, the first and the second time derivatives of the Hubble parameter are given by 


H = 

H = 


Be nt n 2 \ 

(.B + e nt ) 2 ’ 

Be nt (B - e nt ) n 3 A 
(B + e nt ) 3 


(8.59) 

(8.60) 


Furthermore, using in Eq. (2.13) the expression of H derived in Eq. (8.58), we obtain that 
the expression of K 


3se 2nt n 2 X 2 
(B + e nt ) 2 M 2 ' 


(8.61) 


Using in the general expression of Lgo given in Eq. (3.3) the expressions of H and H obtained 
in Eqs. (8.58) and (8.59), we obtain 


Lgo 


(.B + e nt ) 

e nt n 2 X (B/3 + e nt aX) 


(8.62) 


Therefore, we can conclude that the expression of pea with Granda-Oliveros cut-off can be 
written as 


PEA 


3cWA ( BP + e nt aX) 
(B + e nt ) 2 


(8.63) 


Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.58) and (8.63), or equiva¬ 
lently in Eq. (2.22) the expressions of K and pea given in Eqs. (8.61) and (8.63), we obtain 
the following differential equation for F (K), 


dF F 
dK ~ 2K 


a X T /3 


-1 + 


( 3en 2 A 2 A 
\ KM' 1 ) 



eX 


= 0 , 


(8.64) 
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which has a solution given by 


c 2 K ( 2a\ - 2/3 + /3^^ log K) 

F (K) = VKC 5 +-^-r-(8.65) 

SA 

where C 5 represents an integration constant. 

Using in Eq. (2.20) the expression of pea derived in Eq. (8.63) along with the expression 
of H obtained in Eq. (8.58), we can write the pressure pea as 


Pea = - 


c 2 n 2 {R 2 /3 + 3e 2nt aX 2 + Be nt [2aA + /3 (-1 + 3A)]} 
(.B + e nt ) 2 


Therefore, the EoS parameter ujea is given by 


PEA 

uea = - 

Pea 


B(P- Be~ nt (3 - 2aX) 
3A (B/3 + e nt aX) 


At Ricci scale, i.e., for a = 2 and /3 = 1, we obtain 


Lgo 


Pea 


F(K) 

Pea 


uea 


(.B + e nt ) 

y/e nt n 2 X (B + 2e nt A) ’ 

3cWA (B + 2e nt A) 

(B + e nt ) 2 

c 2 K Ux-2 + ^^logK 

VKC 5 + -^^ 

eA 

c 2 n 2 [B 2 + 6e 2nt X 2 + Be nt {7X - 1)] 

(B + e nt ) 2 
B (1 - Be~ nt - 4A) 

_1+ 3A (B + 2e nt X) ' 


( 8 . 66 ) 

(8.67) 

( 8 . 68 ) 

(8.69) 

(8.70) 

(8.71) 

(8.72) 


For a ~ 0.8824 and /3 ~ 0.5016, i.e., for the value of a and /3 corresponding to a non-flat 
Universe, we obtain 


Lgo 


Pea 


F (K) 
Pea 


uea 


(.B + e nt ) 

y/e nt n 2 X (0.50167? + 0.8824e nt A) ’ 

3c 2 e nt n 2 X (0.5016R + 0.8824e nt A) 
(B + e nt f 


(8.73) 

(8.74) 


c 2 K ( 1.7648A - 1.0032 + 0.5016 log K ) 

^KC 5 +-i- - --- l, 

c 2 n 2 {0.5016R 2 + 2.6472e 2nf A 2 + Be nt [1.7648A + 0.5016 (-1 + 3A)]} 

(B + e nt ) 2 

B (0.5016 - 0.5016Re“ nt - 1.7648A) 

_1 + 3A (0.50161? + 0.8824e nt A) ' 


(8.75) 

,(8.76) 

(8.77) 
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For a ~ 0.8502 and /3 ~ 0.4817, i.e., for the value of a and j3 corresponding to a flat Universe, 
we obtain 


Lgo 


Pea 


F(K) 

Pea 


uea 


(B + e nt ) 

y/e nt n 2 A (0.48177? + 0.8502e nt A) ’ 

3cWA (0.4817H + 0.8502e nt A) 
(B + e nt f 


(8.78) 

(8.79) 


c 2 K ( 1.7004A - 0.9634 + 0.4817^^^ log K) 

Vkc, + —i-v-—- L 

e\ 

c 2 n 2 {0.4817H 2 + 2.5506e 2n4 A 2 + Be nt [1.7004A + 0.4817 (-1 + 3A)]} 

(.B + e nt ) 2 

B (0.4817 - 0.48175e“ nt - 1.7004A) 

_1 + 3A (0.48177? + 0.8502e nt A) ' 


(8.80) 

,(8.81) 

(8.82) 


We now consider the Chen-Jing model studied in this paper, i.e., the energy density 
model with the first and the second time derivatives of the Hubble parameter H. Using in 
Eq. (3.7) the expressions of H, H and H obtained in Eqs. (8.58), (8.59) and (8.60), we obtain 
that the expression of pea is given by 


Pea = —— 2 [B [B — e nt ) n 2 a + Be nt n 2 yA + e 2nt n 2 (3\ 2 ] . (8.83) 

Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.58) and (8.83), or equiva¬ 
lently in Eq. (2.22) the expressions of K and pea given in Eqs. (8.61) and (8.83), we obtain 
the following differential equation for F (K), 


dF 

dK 


F c 2 
2 K ~ 3EA 2 X 

< f3\ 2 — (a — qA) 


-1 + 


/ 3en 2 \ 2 
V KM 2 


1/21 


H -(y. 


-1 + 


/3en 2 A 2 
V KM 2 



which solution is given by 


F (K) = + C 6 VK, 

B i 

where Cq represents an integration constant, A\ and B\ are given by 

A 1 = c 2 [— 6 en 2 aA 2 + 2 KM 2 (2a - jA + j3 A 2 ) 

Ien 2 A 2 

+V3KM 2 km2 (-3a + yA) log K , 
Bi = zM 2 A 2 . 


(8.84) 


(8.85) 


( 8 . 86 ) 

(8.87) 
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Using in Eq. (2.20) the expression of pea derived in Eq. (8.83) along with the expression of 
H obtained in Eq. (8.58), we can write the pressure pea as 


where 


Pea = 


A-2 

B 2 


A 2 = c 2 n 2 {—3e 2nt f3\ 3 - B 2 (3a(-1 + A) + 7 A) 

-\-Be nt [o;(—1 + 3A) + A( 7 — 2/3A — 3 7 A)]j, 

B 2 = (B + e nt ) 2 A. 

Therefore, the EoS parameter uiea is given by 

Pea ff{£(3«- 7 A)-e^[« + A (-7 + 2lA)]} 

EA Pea 3A [B 2 a + e 2nt (3\ 2 + Be nt (—a + 7 A)] 

We now consider the scale factor in the intermediate inflation [53, 54]: 

a ( t ) = e Bt8 , 


( 8 . 88 ) 


(8.89) 

(8.90) 


(8.91) 


(8.92) 


where B > 0 and 0 < 6 < 1. With the choice of scale factor given in Eq. (8.92), we obtain 
that the Hubble parameter H, 


H = B6r 1+e . (8.93) 

Moreover, we have that the first and the second time derivatives of the Hubble parameter 
are given by 


H = B(-l + 9)et~ 2+d , (8.94) 

H = B (-2 + 6)(- 1 + 9)8t~ 3+e . (8.95) 

Furthermore, using in Eq. (2.13) the expression of H derived in Eq. (8.93), we obtain that 
the expression of K 


K = 


3 B 2 £t~ 2+26 e 2 

M2 


(8.96) 


Using in the general expression of Lqo given in Eq. (3.3), the expressions of H and H 
obtained in Eqs. (8.93) and (8.94), we obtain 


Lqo 


1 

yjBt- 2 + e 6 [/?(-! + 6) + Bt e a0]' 


(8.97) 


Therefore, we can conclude that the expression of pea with Granda-Oliveros cut-off can be 
written as 


Pea = 3 Bc 2 t~ 2+6 e 


/3(-l + 6) + Bt d aO 


(8.98) 
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Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.92) and (8.98), or equiva¬ 
lently in Eq. (2.22) the expressions of K and pea given in Eqs. (8.96) and (8.98), we obtain 
the following differential equation for F (K) 


dF F 
~dK~2K 




BeQ 


= 0 , 


(8.99) 


which solution is given by 


F{K) 


2 c 2 K 


-e 

( KM 2 \ 2(—1+0) 
\3 B^eff 2 ) 


/3(-l + 6) 2 + Ba6 


Bed 


e 

( KM 2 \ 2(—1+0) 
\3B*e9 2 ) 


+ C 7 \/K,(8.100) 


where C 7 represents an integration constant. Using in Eq. (2.20) the expression of pea derived 
in Eq. (8.98) along with the expression of H obtained in Eq. (8.93), we can write the pressure 
Pea as 


Pea = — 


c 2 {/3(—1 + 6) (-2 + 9 + Wt e 6) + Bt e a6 [-2 + (2 + 3 Bt e ) 0 ]} 

t 2 


( 8 . 101 ) 


Therefore, the EoS parameter ujea is given by 


uea 


Pea _ i t s (—2 + 6) 
Pea 3 B6 


aO 

+ 3/3 ~3Bt e a9 -3fi9' 


At Ricci scale, i.e., for a = 2 and [3 = 1. we obtain 


( 8 . 102 ) 


Lgo 

Pea 


F(K) 

Pea 

uea 


r 

t~ 2+e 9 [{-l + 9) + 2Bt e 9} 
2 t~ 2+e 9 \(-l + 9) + 2Bt e 9 


(8.103) 

(8.104) 


o„2 v ( KM 2 \ 2(-i+0) 


-(-l + fl) 2 + 2R0(j^)^™ 


Be9 


± + VkC 7 ,( 8.105) 


c 2 {(-1 + 9) (-2 + 9 + 3 Bt°9) + 2 Bt e 9 [-2 + (2 + 3 Bt e ) 9] } 


t 2 


t-0(-2 + 9) 
3B9 


29 


3 - 6Bt e 9 - 39' 


(8.106) 

(8.107) 


For a ~ 0.8824 and [3 ~ 0.5016, i.e., for the value of a and [3 corresponding to a non-flat 
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Universe, we obtain 


Lgo 

Pea 


F(K) 


Pea 


uea 


y 7 Bt~ 2 + e 6 [0.5016(—1 + 9) + 0.88245^0] ’ 


Wc 2 t~ 2+e e 


0.5016(—1 + 9) + 0.8824 Bt e 0 


o„2 v ( KM 2 \ 2(-i+0) 
ZC K \3B 2 e8 2 ) 


—0.5016(—1 + 9) 2 + 0.8824 B6 


BeO 




2 

~ x |o. 5016(0 - 1) (0 - 2 + 3 Bt°ej 


+O.8824HU0 


t-%-2 + 0) 
3 B6 


^2 + 3Bt 0 ^ 0 — 2 }, 

0.88240 


1.5048 - 2M72Bt e 6 - 1.50480 ' 


(8.108) 

(8.109) 

0 I 

( KM 2 \ 2(—1+0) 
y 3B 2 e0 2 ) 

( 8 . 110 ) 

( 8 . 111 ) 

( 8 . 112 ) 


Furthermore, for a ~ 0.8502 and /3 ~ 0.4817, i.e., for the value of a and /3 corresponding to 
a flat Universe, we obtain 


Leo 

Pea 


y/Bt~ 2+e 6 [0.4817(—1 + 0) + O.85O2£ e 0] 


3 Bet 


?+-2+e 


0 


0.4817(-1 +0) + O.85O2t 0 0 


(8.113) 

(8.114) 


F(K) 


Pea 

uea 


2 c 2 K 


-9 

( KM 2 \ 2(—1+0) 
\3B' 2 e8 2 ) 


0.4817(-1 + 0) 2 + 0.85020 


0 

( KM 2 2(—1+0) 


Bed 

+Vkc 7 , 

c 2 {0.4817(0 - 1) (0 - 2 + 3 Bt e 6) + 0.8502 1 0 6 [(2 + 3 Bt e ) 0 - 2] } 

t 2 

t~ e (-2 + 6) 0.85020 

W6 + 1.5048 - 2M72Bt e 9 - 1.50480 ' 


(8.115) 

(8.116) 

(8.117) 


We now consider the Chen-Jing model studied in this paper, i.e., the one with the first 
and the second time derivatives of the Hubble parameter H. Using in Eq. (3.7) the expressions 
of H, H and H obtained in Eqs. (8.93), (8.94) and (8.95), we obtain that the expression of 
Pea, 


Pea = 3 c 2 


a (—2 + 0)(—1 + 0 ) 

t 2 


+ Bt~ 2+e -f(-l + 0)6 + B 2 t 


2.-2+28 


P 0 2 


(8.118) 


Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.92) and (8.118), or equiva¬ 
lently in Eq. (2.22) the expressions of K and pea given in Eqs. (8.96) and (8.118), we obtain 
the following differential equation for F (K), 


dF F A 3 
dK ~ 2K ~ 


(8.119) 
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where 


4 -i = 3 


- 2 ‘ 


’(^UV (KM^f 

\B 2 £0 2 ) 


a{-2 + O)(-l + 0) + 3~ 2 e (=h>) B B 


/ KM 2 \ 2 
\B 2 e6 2 ) 


^e) 


ex 


9 -e ( i \<> ( KM 2 \ 2 9 (-i+e) 

7 (-l + 0 ) + 3- 2 ^ *0 (^2 ) 6 


B 3 = B 2 e9 2 . 

The solution of Eq. (8.119) is given by 


F(K) = ^ + C 8 Vk, 

Ba 

where C 8 represents an integration constant, and 

. i e( i ni-er i U 9 / KM 2 \ 4 ^~ 1 + e ) 

A 4 = 2 3 4 v-i+eJ 2 1 -i+6) c 2 K 


\B 2 £0 2 ) 

r t6 a ( 1 \ 29 ( KM 2 \ 2i S ( :r T+e) 

7 ' 726 7 i V30 

2!+3« + 2» (-i,) + 4 (^ 5 ) 

:».w“( r 1 'r 1+2 " 


x 


+ 


+ 


23 1+9 j 7 (- 1 + 6*7 \b 2 e9 2 ) 

7 i \20 7 ) \ 36» 

2 1 + 3 «.B 6 > + 20B [=t+o) 0 + 2B(-^) 0 

3 2l 9 (^+e) B a (2 — 30 + 0 2 ) 


B 2 


l+ 2 1+28 (^) 


—29 


0 2 


B 4 = £. 


( 8 . 120 ) 

( 8 . 121 ) 


( 8 . 122 ) 


(8.123) 

(8.124) 


Using in Eq. (2.20) the expression of pea derived in Eq. (8.118) along with the expression of 
H obtained in Eq. (8.93), we can write the pressure pea as 

PEA = (8.125) 

B 5 

where 


^5 

B 5 


+Bt e 6 

B0. 


ja(—2 + 0)(-l + 0) (-2 + 3 Bt°0} 

7 (-l + 0) (-2 + 0 + 3Bt e 0^j + Bt e p0 ( 


+ 20 + 3 Bt e 0 


}, (8.126) 
(8.127) 
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Therefore, the EoS parameter ujea is given by 


UEA 


PEA _ 2 1 6 

Pea 3 B6 

6('y(-l + e) + 2Bt d p9) 

~ 3{a(—2 + 0)(-l + 0) + Bt e 9 [7 (-1 + 6) + Bt 9 p0}} ' 


(8.128) 


It is possible to analyse matter dominated universe and the accelerated phase of the 
universe using a single formalism. Now for such a model, the Hubble parameter H is given 
by [59, 60] 

H(t) = H 0 + ^, (8.129) 

with Hq and H\ being two constant parameters. From Eq. (8.129), we can easily obtain the 
following expression of the scale factor a(t) 


a(t) = C 9 e Hot t H \ 


(8.130) 


where Cg is an integration constant. Moreover, using Eq. (8.129), we have that the first and 
the second time derivatives of the Hubble parameter, by the following relations: 

H = -|i, (8.131) 

a = (8.132) 

Using in Eq. (2.13) the expression of H given in Eq. (8.129), we obtain the following expression 
for K , 


M 2 


(8.133) 


Using in the general expression of Lqo given in Eq. (3.3), the expressions of H and H 
obtained in Eqs. (8.131) and (8.132), we obtain 

L go = , = — f — . • (8.134) 

yj-pih + a (tHg + fh) 2 

Therefore, we can conclude that the expression of pea with Granda-Oliveros cut-off can be 
written as 


Pea 


3c 2 


-/3Hi + a (tHg + Hi) 2 

T 2 1 


(8.135) 


Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.131) and (8.135), or 
equivalently, in Eq. (2.22) the expressions of K and pea given in Eqs. (8.133) and (8.141), 
we obtain the following differential equation for F (K) 


dF 

F 

3c 2 

' KMP{\eKM*-H 0 ) 2 ' 
3e Hi 

dK 

2 K 


KM 2 


(8.136) 
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which solution is given by 

F (K) = VKC W 

2c 2 [—I8e 2 KM 2 /3 Hq - 27e/3# 0 2 + KM 2 ( e 3 KM 2 p - 9aHi)] 

- 9 TahhI -’ (8 ' 137) 

where Cio represents an integration constant. 

Using in Eq. (2.20) the expression of pea derived in Eq. (8.135) along with the expression 
of H defined in Eq. (8.129), we can write the pressure pea as 

c 2 {—2f3H\ - C tH 0 + #i) [3 t 2 aH 2 + (-2a - 3/3 + 6 taH 0 ) H x + 3 aH 2 ] } . 

PEA = - WWKF) -' (8 ' 138) 

Therefore, we have that the EoS parameter ujea for this case is given by 
ujea = — 1 + 


2taHo 


3 (tHo + H\) 3 \-p Hl + a (t Ho + Hl )‘ 


At Ricci scale, i.e., for a = 2 and (3 = 1, we obtain 
Lgo = 1 


-H^itHo+Hx) 2 

-P- 


-#i + 2 (tH 0 + #i) 


3c 2 


Pea 

F(K) 


-Hi + 2 (tH 0 + H{] 


t 2 

: 2 [-18e 2 AM 2 #o - 27eH% + KM 2 ( e 3 KM 2 - l8Hi)} 
9eM 2 Hi 


(8.139) 

(8.140) 

(8.141) 

(8.142) 

c 2 {-2 Hi - (tH 0 + Hi) [3 t 2 aH 2 + (-7 + 12 tH 0 ) H x + 6 H 2 ] } 

flm + Ho ■ (8 ' 143) 

uea = -1 + —7--- T -,-- (8.144) 

3 (tHo + Hi) 3 ^—Hi + 2 (tHo + Hi) 2 

For a ~ 0.8824 and /3 ~ 0.5016, i.e., for the value of a and (3 corresponding to a non-flat 
universe, we obtain 

(8.145) 


2c 2 [- 

+Vkc 10 , 


Pea 


Lgo 

Pea 

F(K) 

Pea 


^ -0.50i6gi+0.8824(tgQ+gi) J ^-0.5016#! + 0.8824 (tHo + Hi) 

3c 2 [-0.5016tfi + 0.8824 (tH 0 + H x ) 2 


p - ( 8 - 146 ) 

2c 2 [-9.0288e 2 A'M 2 # 0 - 13.5431e#g + KM 2 (0.5016e 3 KM 2 - 7.9416#i)] 

9eM 2 H\ 

+VKC W , (8.147) 

tHtHo + H,) X f - L0032ffl “ {tH " + H ' )X 
[2.6472t 2 #Q + (-3.2696 + 5.2944f# 0 ) #i + 2.6472AT^] } 


uea — — 1 + 


1.7648i# 0 


3 (tH 0 + Hi) 3 _o.5016#i + 0.8824 (tH 0 + Hi 


(8.148) 

(8.149) 
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Furthermore, for a ~ 0.8502 and /3 ~ 0.4817, i.e., for the value of a and /3 corresponding to 
a flat universe, we obtain 


Lgo 


Pea 
F (77) 


Pea 


uea 


-0A817H 1 +0.8502(tH 0 +H 1 ) 2 

-72- 


-0.481777i + 0.8502 (fTT 0 + HiY 


3c 2 


—0.481777i + 0.8502 (t77 0 + Hi) z 


t 2 


(8.150) 


(8.151) 


2c 2 [—8.6706e 2 77M 2 77o - 13.0059^ + KM 2 (0.4817e 3 KM 2 - 7.5618T7i)] 


+VKC 


10 , 


9eM 2 Hi 
x {—0.963477i — (tHo + H\) x 


t 2 (;tHo + Til) 

[2.5506t 2 TT 0 2 + (-3.1455 + 5.1012t77 0 ) 77, + 2.550677?] } 
2 1.7004t77 0 


-1 + 


3 (tHo + H\ 


—0.481777! + 0.8502 ( tH 0 + Hi 


(8.152) 

(8.153) 

(8.154) 


We now consider the Chen-Jing model studied in this paper, i.e., the energy density 
with higher derivatives of the Hubble parameter. Using in Eq. (3.7) the expressions of 77, 77 
and 77 given in Eqs. (8.129), (8.131) and (8.132), we obtain that the expression of pea, 


Pea 


_.tt , 2aH 1 +/3(tH 0 +H 1 ) 3 
A n lE tHo+Hi 


3 £t 2 (77 0 + ^) 2 


(8.155) 


Following the same procedure as in the previous case, we obtain a differential equation for 

F(K), 


F(K) 


Vk 


2c 2 M 2 a 77 3 / 2 c 2 M 2 /3A 3 / 2 

^ 11 H- tt~2 - u - ^ 

e Hi e Hi 

18y/2ec 2 Hq 2 Ha In (77) 


3c 2 \f2MHi)j3 K 

+ 7W1 + 

6\/2ec 2 a 77o 3 In (77) 
Hi 2 M 


77i 2 M 


c 2 M 2 7 77 3 / 2 9c 2 V2MH 0 a 77 

7 

3y/2ec 2 Hq 2 H (3 In (A') 

Hi M 

(8.156) 


where C\\ is an integration constant. Using in Eq. (2.20) the expression of pea derived in 
Eq. (8.155) along with the expression of 77 obtained in Eq. (8.129), we can write the pressure 
Pea as 


Pea = {c 2 [-6a - 2y — 37 2 /377 q + (2/3 + 3y — 6t(3H 0 ) Hi— 

0 ott 2 2t 2 aHl 2taH 0 4a + 27(3a + 7)77 0 

P 1 ” (t77 0 + Hi) 3 ~ (tH 0 + Hi) 2 + 777 0 + 7/, 


x t“ 2 . (8.157) 
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Therefore, we have that the EoS parameter ujea for this case is given by 


loea — {—6o; + 2y + 3t^f3HQ + (—2/3 — 3q + 6f/3-£/o) Hi 

+3BH 2 A- 2t2aH o , 2taH o _ 2[2a + t(3a + 'y)H 0 } 
+ P 1 + {tH 0 + Hi ) 3 {tHv + Hi ) 2 tH Q + Hi 

2aH 1 +(3{tH 0 + H 1 ) 3 Y 


x 


3 -jHi + 


(8.158) 


by 


We can also analyse a q -de Sitter model [62], The scale factor for such a model is given 


a{t) = e q (H 0 t) 


1 + (q- l)H 0 t 


1 

q~ 1 


(8.159) 


With the choice of scale factor, we can derive that the Hubble parameter H along with its 
first and second time derivatives are given, 


H = H 0 [l + H 0 (-l + q)t] 


-i+ 


-1+9 


H = H 2 -1 + 


1 


h = m -2 + 


-1 + q 

i 

-1 + q 


(-! + q) [! + H 0 (-l + q)t\ 2+ - 1+9 , 

1 


-1 + 


-1 + q 


(-! + <?) x 


[1 + H 0 (-l + q)t\ 3+ -i+9. 


(8.160) 

(8.161) 


(8.162) 


Using in Eq. (2.13) the expression of H given in Eq. (8.160), we obtain the following expression 
for K, 


K = 


^ £ Hq [1 + Hq (—1 + q)t] + 1+ 9 
M2 


(8.163) 


Using in the general expression of Lqo given in Eq. (3.3), the expressions of H and H 
obtained in Eqs. (8.160) and (8.161), we obtain 

Lgo = i = = = ^ .(8.164) 

\JH+ H 0 (-l + g)tp 2+3I +l |[1 + H 0 (-l + q)t\~ *+? a - (-2 + <?)/?} 

Therefore, we conclude that the expression of pea with Granda-Oliveros cut-off can be writ¬ 
ten as 

PEA = 3c~H q [1 + Hq{— 1 + q)t] + ~ 1+9 | [1 + Hq(—1 + q)t] ~ 1+q a — (—2 + q)f. (8.165) 

Using in Eq. (2.21) the expressions of H and pea given in Eqs. (8.160) and (8.165), or 
equivalently, in Eq. (2.22) the expressions of K and pea given in Eqs. (8.163) and (8.165), 
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we a differential equation for F ( K ) which solution is given by 


F(K) = 6 c 2 Hl 


K \ 1+ ^7 

ThIJ 


\eHl) 


i 

-1 + 9 


X 


/ 



r / \ —i +9 

\ 

9 I (-2+9) log[X]+log 

K 

+2(— 2-\-q) log 

f K ) 2(—2+g) 
V^O / 

) , 


2(-2+ 9 )2(-l+ 9 ) 4-6(j+2q 2 X 


_ j- T -^4 _ 0 

2(-2+q) 2 (-i+9) (—2 + g) 2 /3 

-1+? 


x M- 2 + \[KC V1 , 


(8.166) 


where C\i is a constant parameter. 

Using in Eq. (2.20) the expression of pea derived in Eq. (8.165) along with the expression 
of H defined in Eq. (8.160), we can write the pressure pea as 


Pea = 


c 2 H% 


[1 + Hq(— 1 + q)t]‘ 


x j-3 [1 + H 0 (—l + q)t] -!+® a 


+(—6 + (7 — 2g)g)/3 + (—2 + q) [1 + FI o(—1 + q)t] 1 +i (2 a + 3/3)j- . (8.167) 


Therefore, the EoS parameter ujea for this case is given by 


ujea = j—3a + [-6 + (7 - 2 q)q] [1 + H 0 (-1 + q)t\ ~'+i /3 
+(—2 + g) [1 + -Ho(—1 + q)t] 1 ~i (2a + 3/3) j- x 
{3a-3(-2 + q)(l + H 0 (-l+q)t)^p} . (8.168) 
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At Ricci scale, i.e., for a = 2 and j3 = 1, we obtain 


Lgo 


Pea 


F(K) 


j-^o [1 + Ho (—1 + q) t] + 1+, x 

{[1 + H 0 (-l + 2 - (-2 + g)}}” 1/2 , 

3c 2 Hq [1 + H 0 (- 1 + q)t]~ 2+ x 
|[1 + 1 + g)t] - 1 +'2 2 — (—2 + g) j , 


6c 2 Hq 


K \ 1+ ^+7 

7hIJ 


\eHfJ 


i 

-1 + 9 


X 


/ 



r / \ —1 + < 7 
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q ( (- 2 + 5 ) log[rf]+log 

K 

[^o\ 

+2( — 2-\-q) log 

/ K \ 2 ( — 2 +9) 

\7ao / 

) 

<? 


2(—2+g) 2 ( — 1+g) 4—6g+2i} 2 x 


(8.169) 

(8.170) 


-l+2g 


/ K \ 2 ( —2+g) 2 ( — l+g) ^ ( — 2 + g)‘ 


x M“ 2 + \/ACii, 


P.E+L = 


= 


c 2 #o 


—9 x / —3 [1 + 77o( — 1 R g)i] 1+9 2 
[1 + 77 0 (—1 + q)t] 2 l 

+(—6 + (7 - 2g)g) + (-2 + g) [1 + if 0 (-l + g)^ ■ 7} , 
{-6 + [-6 + (7 - 2 q)q] [1 + H 0 (-1 + q)t]~^ 

+(— 2 + g) [1 + Hq (—1 + g)f] • 7^ x 
|6 — 3(— 2 + g)(l + Hq( — 1 + q)t) 1-9 | 


(8.171) 


(8.172) 


(8.173) 


For a ~ 0.8824 and /3 ~ 0.5016, i.e., for the value of a and /3 corresponding to a non-flat 































universe, we obtain 
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{[1 + H 0 (-l + q)t]=& 0.8824 - 0.5016(-2 + g)}} ^ , 
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f if A 2(-2+<2) 
V^o7 

) , 


2(-2+g) 2 (—1+9) 4—6q+2q' 2 x 


(8.174) 

(8.175) 


-l + 2q 

K \ 2(-2+q) 2 (-l + q) 


eHl) 


0.8824 - 


(-2 + q) 2 • 0.5016 


M~ 2 + \[KC\i, (8.176) 


c 2h2 


-l+q 

v x (-3 [1 + i? 0 (-l + g)t]=^ 0.8824 
[l + H 0 (-l + q)ty l 

+ [(7 - 2g) g - 6] ■ 0.5016 + (-2 + g) [1 + tf 0 (-l + g)i]^ ■ 3.2696} ,(8.177) 
j—2.6472 + [-6 + (7 - 2g)g] [1 + H 0 (- 1 + g)f] _3 ^ • 0.5016 
+(—2 + q) [1 + H 0 (- 1 + g)*] 1 ^ • 3.2696} x 

|2.6472 -3(-2 + g)(l + 17o(-l + g)i) T ^-0.5016} ' . (8.178) 


Furthermore, for a ~ 0.8502 and /3 ~ 0.4817, i.e., for the value of a and /3 corresponding to 































a flat Universe, we obtain 


Lqo — |t7q [1 “I" 1 + q)t] + H, x 

j[l + 77 0 (-l + g)t]=&5 0.8502 - 0.4817(—2 + g)}} 1/2 

Pea = 3c 2 Hq [1 + 77 0 (-l + g) 7 ] _2+3T +? x 

{[1 + 77 0 (-1 + (/)■£]^ 0.8502 - 0.4817(—2 + g)} , 


(8.179) 

(8.180) 
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[1 + 77 0 (—1 + q)t]“ 

+ [(7 - 2q) q - 6] ■ 0.4817 + (-2 + q) [1 + 77 0 (-l + q)t • 3.1455} ,(8.182) 
| —2.5506 + [-6 + (7 - 2 q)q] [1 + 77 0 (-l + q)t}~^ 0.4817 
+ (—2 + g) [1 + 77 0 (—1 + q)t}^ • 3.1455} X 

j2.5506 -3(-2 + g)(l + TT 0 (-l+ <?)£)^0.4817} ' . (8.183) 


We now consider the Chen-Jing model i.e., the one with the first and the second time 
derivatives of the Hubble parameter 77. Using in Eq. (3.7) the expressions of 77 , 77 and 77 
obtained in Eqs. (8.160), (8.161) and (8.161), we obtain that the expression of pea 

Pea = 3c 2 77q | (—2 + q)(— 3 + 2 q)a + [1 + 77o(— 1 + q)t] ~ 1+IJ x 

(1 + 77 0 (—1 + g)7)^ T +i r /3 - (-2 + q) 7 } x [1 + 77 0 (-l + g)t ]~ 2 . (8.184) 


Following the same procedure as the previous case, we obtain a differential equation for F (77) 
whose solution is given by 


F(K) = \/77Ci3 + 


3c 2 (2 a q 2 


Aa q + 2a — /3 q + f5 + r j) 

- h. 

e 


(8.185) 


where C 13 is a constant of integration. 

Using in Eq. (2.20) the expression of pea derived in Eq. (8.184) along with the expression 
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of H obtained in Eq. (8.160), we can write the pressure pea as 


PEA = -c 2 Hl(l +Hq(-\ +q)t) 2 + i-j x 

2(—2 + q)(— 1 + q )(—3 + 2 q)a + 3 [1 + Hq (—1 + g)t] _1+9 /3 

+(—2 + g)(—3 + 2 q) [1 + H 0 (—1 + g)f] -!+« (3a + 7)— 

(-2 + </)[! + i? 0 (-l + (2/3 + 37)| . 


Thus, the EoS parameter ujea for this case is given by 


loea — — (1 + Hq(— 1 + q)t ) 1 q X 

-J^— 2(—2 + q)(— 1 + g)(—3 + 2g)a + 3(1 + 7/o( — 1 + g)t) _1 +9/3 

+(—2 + g)(-3 + 2g)(l + + g)t)- 1 +9(3a + 7) 

-(-2 + g)(l + l^o(-l + q)t)^ T +^(2f3 + 37) | 

x -^3 (—2 + g)(—3 + 2g)a + [1 + -Ho(—1 T q)t\ - 1 +q x 

(1 + H 0 {-\ + q)t)^+^/3 - (-2 + <7)7 | . 


(8.186) 


(8.187) 


9 Appendix 

In this appendix, we provide the H(z) measurements (in unit [kms” 1 Mpc“ 1 ]) and their 
errors [63]. 
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